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We theoretically study the stability conditions and excitation regimes of hyper-parametric oscillation and Kerr 
frequency comb generation in continuously pumped nonlinear optical resonators possessing anomalous group 
velocity dispersion. We show that both hard and soft excitation regimes are possible in the resonators. Selection 
between the regimes is achieved via change in the parameters of the pumping light. 
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The phenomenon of four-wave mixing (FWM) resuhs in 
hyper-parametric oscillation [1, 2] and frequency comb gen- 
eration [|3H5i] in solid state optical microresonators. The os- 
cillation is closely related to modulation instability lasing 
and parametric oscillation in optical fibers lll0l - [l4ll . 
Four-wave mixing results in generation of frequency side- 
bands in the light leaving the resonator. The sidebands are 
separated approximately by integer numbers of free spectral 
ranges (FSRs) of the resonator and are always equidistant 
{Cjj^ ~ lu = u! — This is because the oscillation re- 

sults from the four photon process hui + hui — >■ hui+ + hCj_, 
where lu+ and tZ)_ are the frequencies of the pump light 
and the generated sidebands, respectively. As the power of 
the pump light is increased multiple sidebands are generated 
and the Kerr frequency comb is produced 

The efficiency of the FWM process depends on group ve- 
locity dispersion (GVD) of the resonator. A resonator with 
nonzero GVD does not have an equidistant spectrum, since 
2a;o — — uj- = cf32Lj-'pgj^/n,Q, where /32 is the GVD (by 
definition, normal dispersion corresponds to f32 > 0), uq is 
the refractive index of the material, assumed to be constant 
in the calculations; c is the speed of light in the vacuum, 
2u!FSR ~ w+ — u!- is the FSR of the resonator; ujo, <^+, 
and are the eigenfrequencies of the consecutive modes. 

In addition to dispersion, the spectrum of the resonator is 
also influenced by self and cross phase modulation effects re- 
sulting from the nonlinearity of the material. In the case of 
anomalous GVD and positive cubic nonlinearity, these effects 
compensate each other so the spectrum of the resonator be- 
comes locally equidistant for a particular power and wave- 
length of light within the resonator The inequidistance of the 
modes cannot be compensated in the case of normal GVD. 
The frequency difference between the modes of the pumped 
resonator influences the FWM and impacts the spectrum and 
threshold of the hyper-parametric oscillation and frequency 
comb generation. 

The Kerr frequency combs excited in optical microres- 
onators are promising for many practical applications since 
their spectral width can span an octave ifTslflql . and their fre- 
quency stability is extremely high [?]. Understanding the fun- 
damental properties of these combs is necessary to achieve 
their optimal performance. This is why, in addition to multi- 
ple experimental investigations, the resonant FWM and KeiT 
comb generation have been studied theoretically. For ex- 
ample, it was shown that additive modulational instability 
ring lasers can operate in both normal and anomalous GVD 



regimes |7]. This result was adjusted to describe the non- 
linear processes in microresonators ifTTUloll . Theories sup- 
porting the idea of preferred generation of the optical Kerr 
comb in resonators possessing anomalous GVD were devel- 
oped lH^im. An analytical expression describing the optical 
pulses generated in resonators with anomalous GVD was de- 
rived 1I22I1 . Finally, generation of Kerr frequency comb was 
investigated via numerical simulations ll23ll24ll . 

The aim of this paper is to analyze the dependence of the 
excitation dynamics of the hyper-parametric oscillation and 
Kerr frequency combs depending on the power and frequency 
of the external continuous wave pump. The dynamics of Kerr 
comb formation was studied previously [24] with approxi- 
mation of small amplitude of the comb frequency harmon- 
ics. This assumption allowed modeling the soft regime of the 
comb excitation. In this work we solved the problem without 
such an approximation, analyze the oscillation onset in the mi- 
croresonators, and s that both hard and soft excitation of the 
oscillation is possible. 

The soft oscillation onset is reached when pump photons 
are not initially present in the resonator. Here, the growth of 
the oscillation sidebands occurs adiabatically. The hard onset 
of the oscillation occurs with a discontinuous jump of the in- 
tensity of oscillation sidebands to a certain finite level, at the 
threshold. The hard excitation occurs only when pump pho- 
tons are initially present in the resonator. The steady state so- 
lution coiTesponding to the hard excitation cannot be reached 
adiabatically. In other words, slow modifications of any pa- 
rameter of the resonator or the pump cannot bring the sys- 
tem to an stable solution requiring hard excitation. Here slow 
means a time shorter than the lifetime of the light confined in 
the resonator. 

We solve the set of nonUnear differential equations describ- 
ing the hyper-parametric oscillation (vf] and frequency comb 
generation 124,1 in steady state, and analyze the stability of 
the solutions in a microresonator possessing a net anomalous 
GVD. While the first order hyper-parametric oscillation (only 
the continuous wave optical pump and the first pair of side- 
bands are involved in the process) is described analytically, 
the higher order oscillations are simulated numerically. 

We found that oscillation sidebands can be much smaller 
than the pump in the case of soft excitation, while they are 
comparable with the pump in the case of hard excitation. The 
sidebands generated in the case of hard excitation produce op- 
tical pulses that travel in the resonator, while the sidebands 
generated in the case of soft excitation are too small to form a 
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pulse. We also find that hard excitation is observed when the 
frequency of the pumping light differs significantly from the 
frequency of the pumped mode, while soft excitation is ob- 
served for the nearly-resonant pumping. We compare results 
of our numerical simulations with the predictions of the an- 
alytical model developed for the comb description ll22ll . and 
find that the simulation predicts formation of wider frequency 
combs and shorter optical pulses relative to the analytically 
found values. The numerical simulations also show that the 
comb spectral width grows slower that the linearized simula- 
tions predict llzsll . 

Let us consider the FWM process that involves only three 
resonant modes. The evolution of the mode amplitudes is de- 
scribed by equations [17] 
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where Fq = i{ujQ — uj)+jo and r± = w±) + 7±, g = 

?iL0Qcn2 / (VtIq) is the coupling parameter iflTll . V is the mode 
volume, 77.2 is the nonlinearity, Fq = {2joP/ {huJo)y^'^ de- 
scribes the amplitude of the continuous wave external pump, 
P is the pump power; A, B+, and B- are the slow ampli- 
tudes of the pump and sidebands respectively. Decay rates 70, 
7_i_, and 7_ reflect both coupling and intrinsic losses of the 
modes. We assume that the modes are overloaded, i.e. the 
loss primarily results from the coupling. 
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FIG. 1: Normalized amplitude of the light within the optically 
pumped mode, and the relative amplitude of the oscillation side- 
bands vs pumping frequency. Stability regions corresponding to soft 
(S) and hard (H) excitation are shown by green solid lines. Unsta- 
ble solutions are depicted by red solid lines. The amplitude of the 
field within the optically pumped mode with no sidebands generated 
{B = 0) is shown by the blue solid line and red dashed line. The 
dashed line stands for the unstable solution for the amplitude of the 
pumped mode when the modulation sidebands are absent. 

To keep a connection between our mathematical models 



and the real physical world we consider oscillations in a mag- 
nesium fluoride whispering gallery mode resonator of 340 /im 
radius (approximately 100 GHz FSR), similar to the resonator 
studied in |26]. The fundamental TE mode of the resonator 
is pumped with 1721 nm light. The modes of the resonator 
have 200 kHz loaded full width at half maximum, and the 
GVD of the modes results in the condition 2ujo — uj+ — uj- = 
—70. The cubic nonlinearity is n2 = 10~^^ cm^/W, and 
the refractive index is no = 1.38. The mode volume is 
V = 1.3 X lO"'^ cm"^ which corresponds to the coupling 
constant g = 1.47 x 10^'^ s^^. By selecting a pump power 
P = 0.235 mW, we obtain (Fo/7o)(g/7o)^^^ = 4for the nor- 
malized pumping constant. It is worth noting that the selected 
value of the pump power is 16 times larger than the threshold 
value needed for the hyper-geometric oscillation to start. 

We solve the set ([T][3]i in steady state, consider a sym- 
metric case, i.e. put 7+ = 7- = 70, so that \B^\ = 
\B^\ — \B\, and present the slow amplitudes of the fields 
as A \A\{1 + SA)exp[i{(j)o + Scjjo)] and B± = \A\{B + 
dB±) cxp[i((/)± + 6cf)±)], where SA, dB+, and SB^ stand for 
the amplitude deviations of the fields related to the drive am- 
plitude, and S(j)o, S(j>^, and S(j>^ stand for the phase deviations 
of the fields. We substitute these expressions into ([T][3l), lin- 
earize the set of equations in the vicinity of the steady state 
solution, and study its eigenvalues. The steady state solution 
is considered as stable if all the eigenvalues are negative. 
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FIG. 2: Normalized amplitude of the light within the optically 
pumped mode and the relative amplitude of the oscillation sidebands 
vs normalized amplitude of the pumping light. The detuning of 
the pumping light from the corresponding optical mode is fixed at 
UJ = ijo — 4.770 (the mode is pumped at its red wind). Dashed blue 
line stands for the undisturbed solution for the light accumulated in 
the pumped mode, B = 0. Red and green solid lines describe unsta- 
ble and stable steady state solutions respectively. 

The results of calculations are shown in Figs. ([T]l and (|2|. 
We have found that there are two regions where the stable 
hyper-parametric oscillation exists. One solution occurs in 
the vicinity of the peak of the resonant curve, shifted due 
to self -phase modulation effect. The oscillation belonging to 
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this branch can be excited adiabatically if one slowly reduces 
the frequency of the pump laser approaching the optical reso- 
nance. It also can be excited if one fixes the laser frequency 
at the red wing of the resonance and then reduces the power 
of the pump (|2]l. The oscillation sidebands are always much 
smaller than the pump amplitude for this solution. The nu- 
merical solution of ([T][3]l shows that the oscillation is excited 
when all the initial conditions are zero. Therefore, we con- 
clude that the stability branch corresponds to the case of soft 
excitation. 

The other stable solution exists when the oscillation side- 
bands have fixed nonzero power. The stability region is local- 
ized in the space of parameters and does not cross the region 
of stable solution for B = 0. Only non-adiabatic change of 
the parameters of the system allow reaching the localized sta- 
bility region. A solution of the initial value problem ([T]|3]l in 
the approximation of zero initial conditions does not reveal 
the stability region. Hence, this stable branch describes the 
oscillation with hard excitation. The localized stable attractor, 
corresponding to the hard excitation regime cannot be discov- 
ered if the original set of equations describing the oscillation 
is solved in the approximation of small sidebands. 
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FIG. 3: Normalized amplitude of the optically pumped mode and 
the relative amplitude of the oscillation sidebands vs pumping fre- 
quency found for the same conditions as used in Fig.[T] The shown 
frequency combs are generated when the parameters of the system 
corresponding to points H and S are selected. 

The analysis that involves only two generated optical side- 
bands is not entirely valid for the description of hyper- 
parametric oscillations in realistic resonators, since a nonlin- 
ear interaction of multiple modes should be taken into ac- 
count. Only the soft-excitation regime that produces small 
oscillation sidebands can be simulated using a limited number 
of interacting modes, as the higher order oscillation sidebands 
are expected to have even smaller power. The hard-excitation 
regime, in which the oscillation sidebands can be more pow- 
erful than the light confined in the pumping mode, cannot be 
strictly understood from the three-mode model. 

To handle this problem we have analyzed the multimode 



regime numerically, deriving a set of equations, similar to ([T]- 
O, for many interacting modes, and solving this set in the 
steady state. With our computing capability generation of op- 
tical frequency combs with mutually interacting 21 optical 
modes can be described. To take the GVD into account we 
assumed that for any pair of sideband modes symmetric with 
respect to the optically pumped mode (ujq — uj^ uj+ — ujq) 
the mode unequidistance is defined by 2wo — uj+ — oj- — 
cl32{u}+ - a;_)2/4no. 
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FIG. 4: Normalized amplitude of the light within the optically 
pumped mode and the relative amplitude of the oscillation side- 
bands vs normalized amplitude of the pumping light. The detun- 
ing is cj = ijjQ — 1570 for plots (a) and (b) (hard excitation), and 
(jj = too + 0.9270 for plots (c) and (d) (soft excitation). The power 
for the first sideband is taken from the carrier. 

The numerical solution revealed excitation regimes simi- 
lar to those found for the case of two optical sidebands (see 
Figs. [3] and lU. A broad frequency comb is generated in the 
occasion of the hard excitation regime. The comb harmon- 
ics have spacings equal to a single FSR of the resonator. The 
frequency comb resulting from soft excitation has harmonics 
separated by three FSRs. There is no single-FSR comb char- 
acterized with soft excitation for the selected pump power. It 
is possible to generate the frequency comb with harmonics 
separated by double or single FSR if smaller optical power is 
selected. 

To demonstrate the hard and soft excitation regimes of the 
oscillation we studied the dynamical behavior of the light in 
the optical modes for zero and nonzero initial conditions (see 
Figs.|5]and|6ll. The outcome of the calculation (the steady state 
solution) is the same for the case of soft excitation (Fig. |5]). 
The oscillation characterized with hard excitation regime oc- 
curs only for nonzero initial conditions (Fig.|6]l. 

Since we can find the phase and the amplitude of the oscilla- 
tion sidebands, we are able to calculate the temporal behavior 
of the field within the resonator We do this for a broad fre- 
quency comb and compare the result with predictions of the 
analytical method described in |22] (Fig.|7]l. As expected, the 
generated harmonics create optical pulses traveling within the 
resonator The numerical modeling and the analytical theory 
predict slightly different behavior for the system. The rea- 
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FIG. 5: Phase diagram for the pump and the first oscillation side- 
band in the case of soft excitation. The blue dashed line corresponds 
to nonzero initial conditions, and the red solid line to zero initial con- 
ditions. The system converges to the same steady state solution with 
different phase of the oscillation sideband. 
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FIG. 6: Phase diagram for the pump and the first oscillation sideband 
in the case of hard excitation. The system does not oscillate if the 
initial field in the pumped mode is zero (red line). The oscillation is 
excited if the initial value of the pump is not zero. 



son is that the analytical model does not take into account the 
finite length of the path of the pulse in the resonator and re- 
quires that the pulse should be generated at the specific detun- 
ing Lj = ojQ — 5.870. The numerical solution is unstable at this 
point, and instead is stable at a neighboring region of detuning 
values (Fig. [3]). 

To conclude, we have developed an analytical model and 
numerical simulations to study the influence of parameters of 
the pump light on the dynamics of KeiT comb generation. We 
found that hyper-parametric oscillation and optical frequency 
comb generation in optical nonlinear resonators can have both 
hard and soft excitation regimes. We present several exam- 
ples of such behavior and show that the hard excitation regime 
leads to formation of short optical pulses in the resonator. 
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FIG. 7: Temporal behavior of the normalized optical field within the 
resonator. Sub-picosecond optical pulses are formed. Inset: the opti- 
cal pulse envelope found from the analytical expression presented in 
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